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Abstract 

We present a general algorithm based on the concept of form- 
invariance which can be used for generating phantom cosmologies. It 
involves linear transformations between the kinetic energy and the 
potential of the scalar field, and transforms solutions of the Einstein- 
Klein-Gordon equations which preserve the weak energy condition into 
others which violate it, while keeping the energy density of the field 
positive. All known solutions representing phantom cosmologies are 
unified by this procedure. Using the general algorithm we obtain 
those solutions and show the relations between them. In addition, 
the scale factors of the product and seed solutions are related by a 
generalization of the well-known a —> a~ l duality. 
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1 Introduction 



Recently obtained astrophysical information ranging from high redshift sur- 
veys of supernovae to Wilkinson Microwave Anisotropy Probe (WMAP) ob- 
servations, indicates that a considerable amount of the total energy of our 
universe might correspond to some dark energy which would be currently 
producing inflation thanks to a large and negative pressure. 

Such exotic fluids may be framed in theories with matter fields that vi- 
olate the weak energy condition pQ. These models were dubbed phantom 
cosmologies, and their study represents a currently active area of research 
in theoretical cosmology. Although phantom cosmologies have been investi- 
gated from different perspectives, here we will only be concerned with par- 
ticular issues related with analytical properties 0. 

At present, the theoretical understanding of the subject is limited, but at 
least one can rely on the motivation for phantom matter which is provided 
by string theory |3J . Precisely, the idea that the origin of dark energy should 
be searched within string theory or another fundamental theory has been 
recently reinforced by the discovery that the holographic principle cannot be 
used to tell whether dark energy is present or not jl]. 

One interesting aspect of phantom matter is that it might make the uni- 
verse end up in a kind of singularity j3] characterized by divergences in the 
scale factor a, the Hubble factor H and its time-derivative H. This possibil- 
ity has not only been shown to exist in general relativity jH], but also in some 
of its generalizations like braneworld cosmology [7] or scalar-tensor theories 
of gravity jH]. 

Another attractive idea put forward recently in connection with phantom 
fields is that they can be used constructing bouncing universes 

In this paper we present an algorithm for generating phantom cosmolo- 
gies with Friedmann- Robertson- Walker (FRW) geometries. Interestingly, the 
method unifies others given in ^1] and [Hj. In fact, here we will prove 
that the triality transformation implemented in [H] which chains inflation- 
ary, cyclic and phantom cosmologies is a form-invariance transformation; 
moreover, such transformation is itself a composition of three simpler form- 
invariance symmetry operations. Thus, that procedure is an application of 
general form-invariance transformations, which have been used successfully 
in ^Tj and for generating new exact solutions. 

Once we discuss the main general features of those transformations we 
concentrate on the linear ones. After that, we show there exist four transfor- 
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mations which preserve the positivity of the energy density, and we demon- 
strate they can be used for generating phantom cosmologies starting from 
standard scalar field universes. We will be mainly be concerned with uni- 
verses represented by scalar field solutions to the Einstein equations, and 
as such the energy density will be sum of a kinetic and a potential term. 
Now, in the usual definition of the energy density of a scalar field the kinetic 
term is taken to be positive, but we will relax this hypothesis and will also 
consider the possibility that it is negative. The transformations we consider 
map between them solutions constructed by taking the two possible signs 
of the kinetic term. This will result in two possible ways of transforming 
those scalar field cosmologies into phantom ones. We will discussed this in 
detail and we will put some accent on comparison with related results which 
appeared earlier in the literature. Finally, we will draw our main conclusions. 

2 The form invariance symmetry 

Our starting point are perfect fluid FRW spacetimes with flat sections given 
by the line element 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). (1) 

Note that, even in such a simple setting as this, many questions remain open 
yet. The Einstein equations are 

3H 2 = p, (2) 

p + 3H(p + p) = 0, (3) 

where p and p are respectively the energy density and the pressure of the 
fluid, and H = a/a. For a different perfect fluid with energy density p and 
pressure p the Einstein equations will take the form 

3H 2 = p, (4) 

p-+3H(p + p) = 0. (5) 

The system of equations ®-© is known to admit form- invariance trans- 
formations |10j . which means there exist symmetry transformations which 
relate it to the system (0)-©- Those transformations are given by 
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p = p{p), 



(6) 
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where p = p(p) is an invertible function. The symmetry transformations 
define a Lie group for any function p(p), which can connect scenarios 
with different evolutions. This will be discussed in the following sections. 

Nevertheless, it is worth making a pair of clarifications before we move 
on. First, form-invariance transformations are defined without imposing any 
restriction on the fluid. Second, they relate solutions to different equations, 
unlike Lie point symmetry transformations, which relate different solutions 
to the same equation. Therefore, form-invariance may be viewed as an un- 
common yet useful equivalence concept. 

The implications of the symmetry transformations will be better under- 
stood if one is able to say how they change relevant physical parameters. We 
will look into this matter from a a general perspective and also from that 
of well-known families of solutions. Under the symmetry transformations 
©-©) the deceleration factor q, which by definition is 



q - -H -, 

a 



(9) 



changes according to 




(10) 



If we consider perfect fluids with equations of state 



P=(7-1)P> P=(7-1)P> (11) 

respectively, then the barotropic index 7 will become 7 and the link between 
them will be 



We have seen that the form-invariance transformations we are concerned 
with require as only input the relation between the energy densities of the two 
fluids. Thus, according to that, we may consider the linear transformation 
between H, p and p generated by 

p = n 2 p, (13) 

where n is a real number. Even though the latter a simple choice, it has got 
much in store. The transformed barotropic index and deceleration factor will 
read 

7=2. (14) 
n 

and 

?=-i + i±i. (is) 

n 

There will be inflation (q < 0) for any n > (q + 1) if n > or for any 
n < (q + 1) if n < 0. The first case was studied in full detail in [TU], whereas 
the second one was partially investigated in jllj . In this case we are going 
to devote ourselves to the second case. 

Finally, for the linear relation (|T3*|) . the transformation rule (JHJ)-© be- 
tween the physical quantities gives 



H \ = \ n \ \ H . (16) 
\ p ) \ n — n 2 n J \ p J 

This means we are actually dealing with a linear transformation acting on 
the variables (p,H,p), which arises because of the particular choice (fT3~j). In 
general other form-invariance transformations will not be linear. Now, the 
second row in Eq. (fT^j) tell us that 

H = nH, (17) 

which leads to 



a = a n . (18) 

The insight gained in this section will be used in the remainder with the 
purpose of generating phantom cosmologies. 
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3 Phantom cosmologies and duality 



When solving the Einstein equations of flat FRW models (0) and Q, one 
looks for solutions of a system of equations in which a does not appear explic- 
itly; H, p and p are the true variables of the problem, and mathematically 
speaking the most natural choice is making assumptions on them, instead of 
the scale factor. 

Here we will consider models such that the Hubble factor is linear in 1/t. 
They correspond to power-law evolutions and have a constant barotropic 
index. Thus, because of the equations of state (fTTj) chosen earlier, we will 
have 

B = W (19) 

H =h m 

where 7 and 7 are the bariotropic index of the two fluids. The scale factor 
of the seed solution associated with (J2U|) reads 

a(t) = (±t) 2/3 \ (21) 

From Eq. (|20|). we have expanding solutions for yt > and contracting 
ones for jt < 0. In order to be able to carry out the following discussion 
for a positive 7, we will take the plus sign for the t > branch of the 
solution and the minus sign for the t < one. The conection between these 
two time regions is made throughout a time reversal symmetry. Since this 
transformation does not change the Einstein equation (0)-©; it is a form- 
invariance transformation. Under that sign choice the two branches will be 
mirror images of each other. Now, replacing Eqs. (|19 |) -([20 |) in (|17|) and 
integrating one arrives at 

^2/37 = | t |2n/37 (22) 

Since in the n < case one has a — > cH"! the transformation may be 
viewed as an extension of the well-known a — ► a -1 duality given in our paper 
[TTj . and which was profusely studied in the pre-big bang scenario [T3]. Now a 
solution of the form of (j21|) with 7 > describes for t > a standard universe 
that expands ever after a big-bang at t — 0, but the t < branch describes 
a contracting solution for which t = represents a big crunch. Now, when 
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Figure 1: The continuous line is a plot of the scale factor as a function 
of time for a standard after-big bang cosmology (t > 0) and its pre-big 
bang continuation. The dashed line represents a new solution obtained by 
applying to the latter a generalized duality transformation a — > a -1 !™'. The 
t < region of the new solution represents a phantom cosmology reaching a 
big-rip at finite time. 



we apply the transformation we obtain a new solution which contracts in its 
t > branch, but expands in the t < one. Precisely, the latter satisfies 
the definition of a phantom universe given in the literature because it is a 
expanding solution that violates the weak energy condition. Moreover, it 
reaches a big rip at finite time, which is a feature inherent to many phantom 
universes but not to all of them. Thus, the phantom universe is dual to the 
universe described by the t < branch of the seed solution. 



4 A unified generation algorithm 

In this section we turn our attention to the scalar field picture; in particular, 
we study how scalar fields transform under the form invariance transforma- 
tions considered above. To this end, we identify the corresponding energy- 
momentum tensor with that of a perfect fluid and write the energy densities 
an pressures in the invariant form 



and 



P(0) = sy + V(0), p(<j>) = s^-V(<j>), (23) 



p{4)=8^ + V($), ptf)=8^-V($). (24) 
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where s and s are two parameters of the scalar field models which can be 
chosen freely. 

Since the energy density and pressure are written as a sum of two terms 
(kinetic and potential energy), there is in principle the possibility of com- 
bining them linearly so that the sums p and p (or p and p) remain constant. 
This generates an additional group structure other than the one introduced 
by the transformation (J13j) . However, since we are only interested in changes 
of sign in the kinetic energy we will restrict ourselves to the discrete symme- 
tries obtained in that case. Therefore, redefining the scalar fields and 0, we 
can select the constant s and s as s = ±1 y s = ±1. Thus, they determine 
the sign of the kinetic terms in the theories described by Eqs. (|23|) and (|24j) 
respectively. 

If we insert (|23|) and (J2"*fj) in Q and © respectively, the original and form- 
invariance transformed Einstein-Klein-Gordon equations will be obtained: 



dV(6) , s 

+ 3H<P + — ^ = 0, (25) 
sdcp 

| + 3i?0 + = 0. (26) 

sdcp 

Applying the form-invariance transformations ©-© and particularizing 
to our choice (|13j) . we find that the kinetic and potential terms transform 
linearly as 

t = -ntf, (27) 
s 

and 

6 2 

V = ns(n - 1)^- + n 2 V. (28) 

Let us see now in more detail how the transformation operates. Under 
the definitions made above we have 

-2H = p + p = s<p 2 . (29) 
and, in consequence, our transformation (fTTjl gives 

p+p = n(p + p). (30) 
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If the seed solution preserves the weak energy condition and considering 
we are in the t < region with n < 0, then the transformed solution will 
violate it and will terminate in a big rip [S] (alternatively see ^H] for big rip 
singularities in non-phantom settings). We concentrate in the remainder on 
the n — — 1 case, which corresponds to p — > p, H — > —if and p+p — > — (/H-p). 
Moreover, the barotropic index of the associate fluid is 

7 = (P + p)/p = -2H/3H 2 = (31) 

P 

and its transformation rule ()12j) reads 7 — > —7. 

Due to the transformation law of the scalar field, the swap between cos- 
mologies that preserve the energy condition and others that violate it can 
occur in two different ways, depending on whether it is the sign of s or 2 
that flips in the transformation 1 . Combining this double possibility with the 
also double choice for the character of the seed solution (s = 1 or s = -1) 
we obtain, from Eqs. (J27|) and (J28j) four possible transformations. We list 
them in Table 1 for the particular case n = —1. 

Table 1: List of n — — 1 transformations. 
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1 




<P 2 + v 


-07p 
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-1 


±0 


<P 2 + v 
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-1 
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±0 


-<P 2 + v 


07p 


-1 


-1 




-<P 2 + v 


-07p 



Now, further insight on how the transformation works can be gained 
from the field equations. We will distinguish between conventional and non 

1 If the sign of <j) 2 is preserved by the transformation the energy density we will behave 
the kinetic energy behaves as a scalar, whereas if that sign changes we will say the kinetic 
energy behaves as a pseudo-scalar. 
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conventional scalar fields by writing the Einstein equations using respectively 
positive or negative kinetic terms. For a conventional scalar field (CSF) the 
set of Einstein equations is 

3H 2 = ^ 2 + V(<j ) ), (32) 

p + 3H<j) 2 = 0, (33) 

whereas for an non conventional scalar field (NCSF) the set of Einstein equa- 
tions is 

ZH 2 = -\tf + V{<P\ (34) 

p - 3H<f) 2 = 0. (35) 

The first transformation in Table I makes 2 — > — 2 , and V — > <fi 2 + V, 
and leaves unaltered the CSF equation set. The second transformation makes 
(f) 2 — > <p 2 and V — > <p 2 + V, so that it turns the set of equations governing 
the NCSF into the set CSF. The third transformation makes <fi 2 — > <fi 2 and 

V — > — 2 + V, and it turns the set of equations governing the CSF into 
the NCSF set. Finally, the fourth transformation makes <j) 2 — > — 2 , and 

V — > — 4> 2 + V, and leaves unaltered the NCSF equation set. 

It is also interesting to look at the implications of the transformations 
from the point of view of the barotropic index 7. A solution to the CSF set 
with a real scalar field will have a positive barotropic index, and its t > 
and t < branches will represent respectively expanding and contracting 
universes with a null scale factor at t — 0. Taking a seed of that sort one 
can construct phantom solutions to the CSF set with an imaginary scalar 
field (using the s — s — 1 transformation) or with a real scalar field (using 
the s = —s = 1 transformation). Specifically, these phantom universes will 
occupy the t < region. In contrast, a solution to the NCSF set having 
a negative barotropic index with a real scalar field and its t > and t < 
branches will represent respectively contracting and expanding universes with 
an infinite scale factor at t — 0, so its t < region branch is a phantom 
universe. Application of the duality transformation gives us a new solution 
such that its t > and t < branches will represent respectively expanding 
and contracting universes with a null scale factor at t — 0. We stress here 
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again that if n 7^ —1 the absolute value of the barotropic index will change 
in the transformation according to 7 = 7/n. 

For illustration we consider now the same families of exact solutions con- 
sidered earlier, that is, the models in which the Hubble factor is proportional 
to 1/1 Choosing all the terms in each equation (j2J) and (f^|) with the same 
degree of homogeneity, the existence of solutions is guaranteed. Hence, the 
corresponding scalar field has a logarithmic dependence on time t and the 
potential depends on the inverse square of t, so that 



0=|ln|t|, 0=|ln|t|, (36) 

and 

V = % V=* (37) 

where A, A, Vq and Vo are parameters related by the transformation rules 
for the field and the potential ()27 )) -(f2"H )) . In addition, they are restricted by 
the respective Einstein-Klein-Gordon equations. Then, inserting Eqs. (|36j) 
and (|S7|) into the Eqs. ([27fl - (j2%]l . we get 

A 2 = —A 2 , (38) 
ns 

9 / 2s T , \ 2ns . . 

Vo = n 2 ( T2 + V ) - — , (39) 

while the composition of (|3^j) with (j3*7j) gives the potentials as functions of 
the fields 

V = V e- A *, V = Voe~ M , (40) 

Requiring that the Hubble factor H, the field <p and the potential V 
given by (fTU|) , and (J3TJ) respectively satisfy the unbarred Einstein-Klein- 
Gordon equations (j2J) and (J25|) we find that 



7= 3l' V °=aAT>- 1 )- (41) 

Finally, replacing the latter in (J39|) we arrive at the complete transformation 
law for the parameter Vq: 
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V> = -g; (Qns - A 2 ) . (42) 

For the seed power law solution a(t) = |t| 2 / 37 or a(t) = \t\ 2s ^ A2 (by using 
Eq. (|4*Tj) ) , the symmetry operations of Table 1 give rise to four different kinds 
of solutions. We list them in Table 2. 

Table 2: Transformations of the power-law solutions in the n = — 1 case. 



s 


s 


4> 




a 




V 


1 


1 






-2/A 2 




> + 


1 


-1 


2,,, 
— In t 




-2/A 2 






-1 


1 


2,,, 
— In \t\ 

A 1 1 


|f 


2/A 2 






-1 


-1 


=F 7 ln|*l 
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The solution generation procedure we discussed in [TT corresponds to 
the first transformation. More specifically, if we take the usual power-law 
cosmologies with a constant barotropic index (the ones used above) and 
apply the transformation s = — 1 = — s we get a family of power-law solutions 
that are, in fact, the late time limit of some generalized phantom cosmologies 
found in ^3] m the context of k-essence cosmologies (see expressions (48) and 
(50) in that reference). If we further set A = A y n = — A 4 /4 in expressions 
()36|) . (|37jl . (|3*5j) . (f59|) and PU|). then we can write those solutions in the form 
in which they were presented in ^3], that is, 

2 A 2 
A = ±-, = ±Aln|t| = ±— 0, (43) 
A 2 

from where it is deduced that 
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A0 = ±^ = \<f), 



(44) 



and finally 

K, = y(3A 2 + 2), (45) 

V = V e- A * = V e^ /X = V e~ x *. (46) 

Then, the solution presented in ^3] corresponds to a particular appli- 
cation of the second transformation type. Actually, in his paper ^3], the 
author only presented models with potentials that decrease as <fi increases, 
which follow from taking the lower sign in (|43|) and (|46|). 

Interestingly, using the concept form-invariance, and with the help of 
the parameters syswe have unified the procedures already existing in the 
literature, and as a byproduct we have obtained the alternative procedures 
associated with the third and fourth transformations. 



5 Conclusions 

Our main goal here has been giving a general prescription for generating 
phantom cosmologies which unifies those obtained from the Einstein-Klein- 
Gordon with both sign of the kinetic term. The algorithm exploits form- 
invariance transformations, which have been discussed earlier in relation to 
cosmology as part of a long-term project. After some revision of the most 
important properties of the transformations, we have concentrated on linear 
transformations that multiply the Hubble factor by a negative constant factor 
n, so that contracting solutions will be transformed into expanding ones and 
vice versa. The transformation in the scale factor will be a generalized version 
of the a — > a^ 1 duality investigated in [TT] and appearing in the pre-big bang 
scenario, specifically our transformation will consist in a — > a~' ra '. 

Let us assume n — — 1 and that the seed solution is the typical contracting 
(expanding) solution in the t < (t > 0) region ending (beginning) in a 
singularity at t — (see Fig. 1). Given that a — > a^ 1 , the t < branch of 
the transformed solution will have a expanding character and it will reach a 
big rip at t — 0, thus representing a phantom cosmology. The t > branches 
of the transformed and seed solutions are the standard duals of each other. 
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For a constant n, the form-invariance transformations we consider can be 
viewed as group of uniparametric linear transformations acting on the sets 
of variables (p, H, p) or (0 2 , H, V), where 2 and V are the kinetic term and 
the potential of a conventional scalar field. Although we fix the value of 
the parameter of the group n, some discrete symmetry remains associated 
with the possibility of choosing a positive or a negative sign in kinetic term 
in the effective energy density of the scalar field. In addition, two different 
symmetry transformations can be constructed between the CSF equation set 
and the NCSF equation set. Combining those possibilities, we obtain four 
transformations that revert the sign of the barotropic index, and two of them 
generate phantom cosmologies. 

In conclusion, the procedures to transform conventional scalar field cos- 
mologies into phantom ones which existed up to now in the literature can be 
viewed as an application of form-invariance transformations. 
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